Anti-atherosclerotic effects of naringenin and quercetin from Folium Artemisiae argyi by attenuating Interleukin-1 beta (IL-1β)/ matrix metalloproteinase 9 (MMP9): network pharmacology-based analysis and validation

Effective components and related target genes of Folium Artemisiae argyi were screened from Traditional Chinese Medicines for Systems Pharmacology Database and Analysis Platform. The therapeutic targets of atherosclerosis were searched in the MalaCards and OMIM databases. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed in WebGestalt online and verified according to ClueGo and Pedia apps in Cytoscape. Then, the protein-protein interaction network was analyzed using the STRING database and constructed using Cytoscape. Differential expression of target genes was identified in GSE9128 and GSE71226 by GEO2R. And then, molecular docking was performed using the Molecular Operating Environment. Finally, we validated the protein expression of Interleukin-6 (IL-6)/IL-1β /MMP9 by qRT-PCR and Western blot in Raw264.7 which was induced by LPS. A total of 232 potential target genes and 8 ingredients of Folium Artemisiae argyi were identified. Quercetin and naringenin are potential candidate bioactive agents in treating atherosclerosis. Vascular endothelial growth factor (VEGFA), MMP9 and IL-1β could be potential target genes. KEGG analysis demonstrated that the fluid shear stress and atherosclerosis pathway play a crucial role in the anti-atherosclerosis effect of Folium Artemisiae argyi. Gene Expression Omnibus (GEO) validation demonstrated that VEGFA was downregulated, while MMP9 and IL-1β were upregulated in patients with atherosclerosis. Molecular docking suggested that only MMP9 had a good combination with quercetin. The cell experiment results suggested that naringenin and quercetin have strong anti-inflammation effects, and significantly inhibit the expression of MMP9. Practical Applications Artemisiae argyi is a traditional Chinese herbal medicine that has been widely used for its antibacterial and anti-inflammatory effects. This research demonstrated the bioactive ingredients, potential targets, and molecular mechanism of Folium Artemisiae argyi in treating atherosclerosis. It also suggests a reliable approach in investigating the therapeutic effect of traditional Chinese herbal medicine in treating Atherosclerotic cardiovascular disease (ASCVD). Supplementary Information The online version contains supplementary material available at 10.1186/s12906-023-04223-1.

Anti-atherosclerotic effects of naringenin and quercetin from Folium Artemisiae argyi by attenuating Interleukin-1 beta (IL-1β)/ matrix metalloproteinase 9 (MMP9): network pharmacology-based analysis and validation Background Atherosclerotic cardiovascular disease is a common pathology present in many cardiovascular diseases, arising from the obstruction of coronary vessels due to atherosclerosis or thrombosis.Cardiovascular diseases (CVDs) were the leading cause of death in non-communicable diseases according to the World Health Organization over the past decades [1].CVD and stroke have become the two greatest causes of burden of disease in high-income countries.The risk factors of CVD vary and include sex, smoking, alcohol intake, and deficiencies in social relationships [2].The causes of utmost concern in the general population are dyslipidemia and inflammation.Thus, statins and other drugs are widely used to reduce lipid levels and inflammation in hypercholesterolemia and other CVDs [3].Although the effect is obvious, the side effects in some patients are also unavoidable [4].The use of Traditional Chinese Medicine combination might be a potential supplementary treatment.
Artemisiae argyi is a Chinese herbal medicine containing many bioactive compounds, such as flavonoids, glycosides, sterols, and essential oils [5].It has been widely used for the treatment of infections, cancers, and other inflammatory diseases [6].Inflammation is also a major cause of CVDs [7].Therefore, we explored the possible treatment of ASCVD using Artemisia argyi according to a pharmacology-based network analysis method [8].
Network pharmacology [9] is a novel paradigm that integrates the concepts of network science and pharmacology, provides numerous valuable advantages in the research of drug discovery and development.The traditional medicine pharmacology network prediction analysis is a method involving the pharmacogenomics and therapeutic mechanism of traditional Chinese medicinal herbs and/or formulae and the potential target genes and/or drugs [10].The comprehensive investigation of the relationships among drugs, target genes, and diseases are possible because of the rapid development of bioinformatics and pharmacology [11].

Pharmacokinetic absorption, distribution, metabolism, and excretion (ADME) screen
The ADME criteria of Folium Artemisiae argyi were extracted from the TCMSP database.Drug-likeness (DL) and oral bioavailability (OB) were selected to identify the bioactive ingredients of Folium Artemisiae argyi.OB is the percentage of an oral dose capable of producing pharmacological activity [14].DL is an indicator for determining the similarity or likeness of a compound that can help in determining whether a compound has a therapeutic effect or not [15].

Targets of compounds searching
Information on the compounded ingredient target genes was obtained from the TCMSP database, and the Drug Bank [16] (https://go.drugbank.com/)database was also used for determining the comprehensive drug targets of all ingredients.The related target genes of atherosclerosis were searched from the Mala Cards (https://www.malacards.org/) and OMIM (https://omim.org/)databases.The target genes of compounds were collected according to the Similarity ensemble approach (SEA) [17] online database (http://sea.bkslab.org/).

Protein-protein interaction (PPI) network
The overlapping genes of AS and the compounds were selected as the hub genes and analyzed using the database STRING [18] (https://string-db.org),which could provide the PPI network results.The Cytoscape [19] (https:// cytoscape.org/)software is widely applied to pharmacology studies in network construct and visualization.

KEGG analysis and enrichment
KEGG database was established by the Kanehisa Laboratory in 1995 and is typically used in pathway analysis and annotation in network pharmacology.We used WebGestalt [20] (WEB-based Gene Set Analysis Toolkit, http://www.webgestalt.org/)for KEGG pathway analysis, which is a functional enrichment analysis web tool.Then, the interactions between genes and pathways were validated by ClueGo and Pedia apps in Cytoscape.

GEO validation
Candidate target genes were identified in the GEO database (GSE9128, GSE71226), including total RNA expression data from coronary heart disease in human.GEO2R was used to identify the differentially expressed genes (DEGs), p ≤ 0.05, and |log FC| > 1 were the screening limitations [21].

Molecular docking
Molecular docking is a crucial technology in proteins and small compounds.It is performed using the Molecular Operating Environment (MOE, v2019.0102)software to validate interactions between compounds and target proteins.The 3D structure of targets was obtained from the Protein Data Bank (PDB, http://www.rcsb.org)and imported into MOE to perform molecular docking after protein structure preparation.The structure of participant compounds was obtained from PubChem (https:// pubchem.ncbi.nlm.nih.gov).

Compounding ingredients of Folium Artemisiae argyi
We input "Folium Artemisiae Argyi" as an "herb name" to search the ingredients of the compound.A total of 135 items were obtained, and only 9 ingredients were included after screening by OB ≥ 30% and DL ≥ 0.18 in this study (Table 2).The target genes of the nine ingredients were collected from the SEA (Similarity ensemble approach, http://sea.bkslab.org/)database, which is a database that can be searched for chemical formulas according to their ingredients.After selecting genes from humans and eliminating the duplicate genes, 8 ingredients and 232 genes were included.The network of compounds to target genes was constructed using Cytoscape (Fig. 1).

Atherosclerosis-related target genes
We searched the Mala Cards database online using the keyword "atherosclerosis" and 81 target genes were selected.We searched the Online Mendelian Inheritance in Man (OMIM) database with the same keywords, and 269 genes were collected.The overlapping genes among "atherosclerosis-related genes" from Mala genes and OMIM database and the "compound ingredients target genes"; finally, eight hub genes of the nodes with the highest connections that appear in the network were identified (Fig. 2).
Fig. 1 Protein-protein interaction network of common targets of Folium Artemisiae argyi and atherosclerosis.Red circle: Folium Artemisiae argyi; Purple hexagon: ingredients of Folium Artemisiae argyi after screening by OB ≥ 30% and DL ≥ 0.18; green diamonds: target genes (human) of the ingredients; blue diamonds: eight hub genes overlapped among "atherosclerosis-related genes" from the Malagenes and OMIM databases and "compounding ingredients target genes."

Network construction of protein-protein interaction (PPI)
The eight hub genes were input into the online tool "STRING", and the PPI network was constructed by the limitation: "minimum required interaction score, (confidence = 0.500); max number of interactors, (1st shell ≤ 20 interactors, 2nd shell ≤ 20 interactors)"; in total, 48 interactors were collected.Visualization of common targets and related signaling pathways was created, larger node in the graph means higher degree value and greater likelihood to the targets (Fig. 3).KEGG pathway analysis of the 48 genes was performed using the online web tools WebGestalt.The top 20 pathways are listed in Table 3; Fig. 4. The results were validated using ClueGo + Pedia apps, and the "Fluid shear stress and atherosclerosis" pathway (Fig. 3), including the three genes (IL-1β, MMP9, VEGFA), were the target genes in the eight hub genes.

GEO validation
Analysis of GSE71226 and GSE9128 expression data of atherosclerosis group revealed that VEGFA was downregulated, while MMP9 and IL-1β were upregulated (Table 4, Supplementary Tables 1 and 2).Thus, the three genes might be the candidate therapeutic targets of Folium Artemisiae argyi in the clinical treatment of atherosclerosis.

MOE docking
Molecular docking was performed to validate the interaction of the target protein IL-1β (PDB code: 2nvh), MMP9 (PDB code: 2ow1), and VEGFA (PDB code: 1mkk) and the related participant compounds (MOL000098, and MOL001040).Visualization of docking results of receptors and ligands were presented in Figs. 5 and 6.The results of Fig. 5A-C showed that quercetin could interact with MMP9, IL-1β, VEGFA. Figure 5A shows that quercetin can interact with Leu132 and Glu130 in MMP9, Fig. 5B shows quercetin interact with Lys97 in IL-1β, and Fig. 5C shows quercetin interact with Asp34, Leu32 and Glu30 in VEGFA.Figure 6 A represent naringenin interact with His401 and Leu418 in MMP9, Fig. 6B shows naringenin interact with Lys63, Lys65 and Gln38 in IL-1β, Fig. 6C shows naringenin interact with Ser50 and Phe47 in VEGFA.The binding energy score is shown in Table 5. Lower binding energy means higher affinity between the receptor and the ligand, and the conformation is more stable.Generally, binding energy which is less than − 5 kcal/mol indicates a good binding activity [22].The molecular docking results predicted that quercetin and naringenin from Folium Artemisiae argyi could be effective in atherosclerosis therapy by targeting MMP9.

Quercetin and naringenin suppressed LPS-induced proinflammatory cytokines
To investigate the effects of quercetin and naringenin on anti-inflammation, Raw264.7 were stimulated with LPS in the presence or absence of quercetin (10, 20, 50 µM) and naringenin (10, 20, 50 µM) for 24 h.As shown in Figs.7 and 8, the mRNA and protein expression of IL-6 and IL-1β were significantly increased with the LPS treatment (P < 0.0001), which were inhibited by quercetin (10, 20, 50 µM) and naringenin (10, 20, 50 µM).These results provided evidence that quercetin Fig. 2 The Venn diagram revealed the overlapping genes among "atherosclerosis-related genes" from Mala genes and OMIM database and the "compound ingredients target genes".Gray circle: "atherosclerosis-related genes" from the Mala genes database; orange circle: "atherosclerosis-related genes" from the OMIM database; and blue circle: "compounding ingredients target genes" from the SEA database and naringenin have a strong inhibitory effect on proinflammatory cytokines.

MMP9 might be a therapeutic target of Folium Artemisiae argyi in ASCVD
The mRNA and protein expression level of MMP9 had a significant increase after treatment with LPS in the Raw264.7,quercetin and naringenin could significantly decrease the mRNA and protein expression level of MMP9 (Fig. 7 C, D and Fig. 8), suggesting that lowering the expression of MMP9 might be the therapeutic effect of quercetin and naringenin in the treatment of ASCVD.

Discussion
Multiple factors are associated with cardiovascular diseases [23]; a high low-density lipoprotein cholesterol (LDL-C) concentration of plasma and inflammation are the major factors that cause atherosclerosis.Lowing cholesterol uptake and increasing cholesterol efflux could attenuating the inflammation of atherosclerosis plaque and even promote plaque regression [24].Although statin treatment in lowering LDL-C has achieved a relatively optimistic result, its benefits are limited by adverse effects to the liver and others, and further effective drugs for treating ASCVD should be sought.Inflammation plays an important role in theprocess of atherosclerosis, the Canakinumab Anti-inflammatory Thrombosis Outcomes Study (CANTOS) showed that targeting inflammation in atherosclerosis in clinic was effective [25].
In the present work, we constructed a network of bioactive compounds and the molecular targets of Folium Artemisiae argyi that overlapped the target genes between atherosclerosis and related ingredients of Folium Artemisiae argyi.Finally, eight hub genes were identified, and IL-1β, VEGFA, and MMP9 genes in the fluid shear stress and atherosclerosis pathway are the most likely target genes in treating atherosclerosis.The results of the GEO database (GSE71226,   GSE9128) validation revealed that IL-1β and MMP9 expression was upregulated, and VEGFA was downregulated significantly compared with controls.The expression values of VEGFA from GSE9128 in the control group were higher compared to the Ischemic cardiomyopathy (ICM) group, coincident with early research that increased the expression of VEGFA might be a potential therapeutic method for ICM [26].MMP9 involved in the matrix-metalloproteinases family has been implicated in regulating matrikines.Given their ability to alter cellular migration and mitogenesis, matrikines have been implicated in inflammation, wound repair, and atherosclerosis [27].MMP9 plays a role in inflammation and is upregulated in a lipopolysaccharide (LPS) model of corneal inflammation [28].And the molecular docking results also suggested that MMP9 has better interactivity with quercetin, the experiment results in Raw264.7 were also providing evidence that quercetin and naringenin could decrease the expression of MMP9 and suppressed the expression of pro-inflammation cytokines IL-6 and IL-1Β.
In our study, the LPS induced inflammation in Raw264.7 also elevated the mRNA and protein expression level of MMP9, as while treated with quercetin (10, 20, 50 µM) and naringenin (10, 20, 50 µM) could significantly decrease its expression (Fig. 7 C, D and Fig. 8).And these results suggested quercetin might have the effect of steady atherosclerotic plaque stability by inhibiting MMP9 expression.
IL-1β is a member of the IL-1 family cytokines; it is an immunomodulatory signaling molecule and thus acts as a central mediator [29].The Canakinumab Anti-Inflammatory Thrombosis Outcome Study trial also provided proof for the inflammation hypothesis of atherosclerosis [30], and IL-1β inhibition highlighted the potential of anti-inflammatory therapies to improve the clinical outcomes of CVDs.The results of our research also presented the suppression of quercetin and naringenin to pro-inflammation cytokines in IL-1β and IL6 (Figs. 7A  and B and 8).
The MMP9 and IL-1β-related major ingredients of Folium Artemisiae argyi were quercetin and naringenin.Quercetin, one of the ingredients of Folium Artemisiae argyi, has shown a wide range of biological actions in anti-inflammatory and antiviral activities in vitro and in some animal models [31].The ability of inflammation to promote atherosclerosis has been elucidated in molecular and cellular pathways by numerous experimental works [32].Quercetin is a kind of flavonoid, and a prominent dietary antioxidant present in fruits, vegetables, and herbal medicines, it plays a role in attenuating atherosclerosis by alleviating inflammation and improving nitric oxide (NO) bioavailability [33].
Naringenin is also one of the natural flavanones in Folium Artemisiae argyi, and animal models have demonstrated its therapeutic potentials in treating inflammation-related diseases, such as atherosclerosis [34].Naringenin suppresses inflammatory cytokine production during transcription and post-transcription; it not only inhibits cytokine mRNA expression but also promotes lysosome-dependent cytokine protein degradation [35].Thus, the anti-atherosclerotic activity of naringenin is due to its high anti-inflammatory effects [36].
Inflammation is an important driver of atherosclerosis, and cellular inflammatory changes actively contribute to atherosclerosis progression [37,38].The therapeutic effect of the inflammatory pathway targets helped improve the outcomes of patients with cardiovascular diseases.The anti-inflammatory effects of the ingredients from Folium Artemisiae argyi were obvious.Consequently, Folium Artemisiae argyi has potential beneficial effects in atherosclerosis therapy through its anti-inflammatory activities.However, our research has limitations in investigating the mechanism of Folium Artemisiae argyi used in treating atherosclerosis.And its application to clinical medicine in the future should be determined through extensive experiments in vivo and in vitro.

Conclusions
In the present study, we performed network pharmacology-based prediction, molecular docking, and GEO database validation to verify the potential targets of Folium Artemisiae argyi through related bioactive ingredients in treating atherosclerosis.And the validation in the LPS-induced inflammation model of Raw264.7 also offered evidence that quercetin and naringenin have the anti-inflammation effect and suppressed the expression of MMP9.We demonstrated that the anti-inflammatory and keeping the atherosclerotic plaque stable ability of Folium Artemisiae argyi may be the main direction in atherosclerosis therapy in the future, which also provided a practicable application for the analysis of traditional Chinese medicine in the clinical treatment of diseases.

Fig. 3
Fig.3Network of the main pathways and targets of Folium Artemisiae argyi treatment of atherosclerosis."Fluid shear stress and atherosclerosis" pathway[39] and related genes. in the center of the circle were the genes related "Fluid shear stress and atherosclerosis" (linked with red lines)

Fig. 4 Figs. 6 Figs. 5
Fig.4 The KEGG pathway enrichment analysis results.Top 20 pathways analyzed in WebGestalt, arranged by the enrichment ratio from high to low

Table 4
GEO validation using GEO2R

Table 5
The binding energy score of ingredients with MMP9, IL-1β and VEGFA s (Kcal/Mol).